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1 Introduction 



Dynamical fermion simulations require the calculation of statistical averages with 
the pure gauge measure modified by the fermion determinant. Its evaluation is at 
present the most time consuming step of lattice QCD simulations. In spite of the 
effort spent to evaluate quantum averages in the full theory, the results from these 
calculations when compared with those of the quenched approximation, where the 
fermion determinant is set to unity, appear to amount mainly to a shift in the values 
of the bare coupling (3. The size of the effect can be estimated by requiring that the 
lattice spacing of the pure gauge theory (G) matches the one of the full theory (GF 
= Gauge + Fermions) in perturbation theory: 

a GF {(3 GF ) = a G {(3 G ) 

a { A = e~ 12b o (6bi/(3)~^ i = GF, G (1) 

where A is the same in the two theories because physical quantities are set to com- 
mon values and b and b\ are the coefficients of one and two loops beta function. 
What remains of fermion dynamics after the shift appears to be a perturbative cor- 
rection to the quenched result: in this paper we argue that such a correction can be 
extrapolated from a theory with negative flavour numbers where the fermion fields 
are replaced by fields with Bose statistics and fermion action that we call bermions. 
The extrapolation from negative flavour numbers was proposed long ago in ref. [1] 
for two dimensional models, but the method worked only for heavy quark masses. 
The results described in this paper have been obtained with a Wilson fermion ac- 
tion corresponding to even flavour numbers. This restriction is due to the necessity 
of dealing with a positive definite quadratic operator for the bermions which are 
thermalized together with gauge fields in the Monte Carlo simulation. It is the lat- 
ter possibility which makes the method very competitive with respect to the best 
fermion determinant algorithms: we have estimated, for the choice of volume and 
parameters of this paper, a gain by a factor approximately one hundred with respect 
to a hybrid Monte Carlo simulation. 

In the first part of the paper after a presentation of some technical details we 
discuss how the extrapolation works when the shift is expected to be perturbative, 
i.e. for heavy quark masses or in the asymptotic scaling region; in the second part 
we show how, even for a non perturbative shift, the extrapolation remains smooth 
when it is performed at constant values of the lattice spacing and of the renormalized 
quark mass. 
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2 The bermion action 



Following Ltischer in ref. [2] we use the action 

s[u,4>] = s g [u] + Y.<t>\ x )Q 2( t>{ x ) (2) 

X 

where S g is the standard Wilson action for the gauge part, Q is an operator related 
to the lattice Dirac operator D by 

Q = 75 {D + m)/M, (3) 

and (j)(x) is the bermion field. 

The constant M is related to m by: 

M = 8 + m (4) 
Introducing the Wilson hopping parameter 

K = (8 + 2m)" 1 (5) 

and the constant 

c = {2KM)- 1 = (1 + 8K)- 1 (6) 

one obtains the result: 

[Q<f>](x) = Co l5 {</>(x) - K " l^( x + A*) 

jii=0 

+Ul(x- f i)(l+ 7„M*-M)} (7) 

The form of the bermion action is suitable for a simulation with a standard heat- 
bath algorithm. Actually we did use a hybrid algorithm which involves a mixture 
of heat -bath and of overrelaxation on the gauge fields and on the bermion fields. In 
the SU(3) case we used the Cabibbo-Marinari pseudo-heat-bath algorithm. A full 
sweep of the lattice is composed of one heat -bath and iVj, overrelaxation sweeps over 
the bermion fields, followed by one pseudo-heat -bath and N g pseudo-overrelaxation 
sweeps over the gauge fields. The parameters iVj, and N g are optimized in order to 
keep auto-correlation times as low as possible. For hadron spectroscopy, we need 
the inversion of the quark propagator that we perform with the conjugate gradient 
algorithm every 1000 full sweeps. 
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3 Perturbative dynamical fermions 



When the shift is perturbative, the connection between negative and positive flavour 
numbers can be obtained by a fermion loop expansion [3] and can be expressed in a 
simple analytic form. If we represent the fermion average as: 

F _ J D[A]e- s WQ[A]e(^W 
[ )n f J D[A]e- s W n tM*W ^> 

the corresponding bermion average is expressed by: 

J D[A]e- s WQ[A]e-^W 

\ U /n b j D ^ e -S[A] e -n b A[A] ^ > 

where rtj and n& are both positive and integer numbers and A = Tr(\og(Q 2 )). From 
bermion averages one can produce fermion averages which are correct up to a definite 
number of fermion loops. For example, an average valid to one loop can be written 
as: 



(on ~ r ^ do) 

Similarly, one can build expressions corrected up to a higher loop number: 



(O) 



2 loop (((9) ) 3 ((9)f 
2 ~ ((O)ff 



F 



F 3 1oo P (((9) ) 4 (((9)f) 4 



The equations above are expected to converge rapidly to the unquenched result 
for heavy quark masses. In table 1 we give the value of the average plaquette for the 
SU(2) case at a value of (3 far from the asymptotic scaling region at two values of the 
Wilson parameter K. For the first rather low value, i.e. for heavy quark masses, the 
jump in the plaquette is obtained correctly by using the simple one loop formula in 
eq.(10). For the second K value, i.e. for lighter quarks, the jump increases and the 
one loop approximation becomes inadequate. The inclusion of higher loops lowers 
the discrepancy with the right value but the error of the estimate increases. This 
loss of precision can be understood by considering that equations (11) imply a ratio 
of several correlations with a natural increase of the resulting statistical fluctuations. 

The steeper flavour number dependence for lighter quarks can be seen from 
figure 1 where we plot the average plaquette for (3 = 1.75 and for different values 
of the quark mass. The effect can be attributed to the increase in the shift of the 
(3 values when quarks are lighter and make the contribution of fermion loops fully 
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active; the asymmetry with respect to zero flavour values reflects the different slope 
of the average plaquette at (3 values lower and higher than the one of the pure gauge 
case. 

A fast convergence is also expected in the asymptotic scaling limit at very 
large values of (3. We report the perturbative expressions in SU(3) for the average 
plaquette (P) and for the quark mass: 

, 4 2 1.2208 n f 

1 " < p > = p + — " °- 0588 ^ 

1 1 N , , 1 5.457-3, , , 

-, = ln(l + S - il -) = ln(l + --4 + ^^) (12) 

In eq.(12) the fermion mass in plaquette expansion is taken to be zero (see ref. [4]). 
From now on our results refer to the SU(3) case only. 

We first explore the case of (3 ~ 24 on a 8 3 • 16 lattice, deep in the asymptotic 
scaling region. In figure 2 we report the lxl and 2x2 plaquettes as a function 
of the flavour number: the result are mildly dependent upon the flavour number. 
Figure 3 is a blow up of the behaviour of the 2x2 plaquette: in such an enlarged 
scale the flavour dependence is visible and appears well fitted by a straight line. 
The lowest order of the perturbative expansion proportional to the flavour number 
is adequate at such a high (3 value and provides a nice test of the bermion Monte 
Carlo. 

At lower values of (3 the perturbative expansion breaks down: Parisi has sug- 
gested [5] that its validity may be extended by a suitable choice of the expansion 
parameter eliminating tadpole type corrections which are typical of lattice perturba- 
tive calculations and responsible for the large ratio of the lattice A parameter with 
respect to, for example, the one of dimensional regularization. The original proposal 
was to adopt an expansion parameter related to the value of the lxl plaquette. 
This suggestion was later refined by Lepage and Mackenzie [6] who proposed an 
improved definition of the coupling constant still related to the lxl plaquette. In 
the same paper they proposed an expression for the quark mass in terms of the im- 
proved coupling evaluated at a different and optimized scale. Comparing different 
theories at the same value of the improved coupling requires a readjustement of the 
corresponding (3 values to obtain a matching of the lxl plaquette from which the 
improved coupling is extracted. 

The strategy of matching a physical quantity like the improved coupling at the 
same lattice spacing scale can be further generalized into the recipe of comparing 
and then extrapolating theories with different flavour content at fixed values of the 
renormalized quantities. This procedure reabsorbs the unphysical shift in the bare 
couplings and allows to extrapolate only the genuine physical effects. In QCD with 
flavour degeneracy there are two independent quantities which need a renormaliza- 
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tion: the coupling constant and the quark mass. We then need two independent 
quantities to match in order to keep renormalized quantities fixed. At perturbative 
level these can be identified with the improved coupling and the quark mass defined 
in ref. [6]. The effect of the matching is visible also in the perturbative region: in 
figure 4 we report the variation with the flavour number of the 2x2 plaquette to be 
compared with the one in figure 3. The behaviour is still well described by a linear 
fit, but the variation itself has been reduced by a factor 4. 



4 Non perturbative dynamical fermions 

The procedure described in the last section can be fully appreciated at moderate 
values of (3: we present the results at fixed bare parameters for (3 = 5.4 and K = .162 
on a 16 3 • 32 volume in figure 5, where now the flavour dependence is visible also 
with a poor scale resolution. After the matching described above, the same figure 
turns into figure 6 where the variations are much smaller and linear. A magnified 
scale allows to better appreciate the change in the behaviour: we give in figure 7 the 
2x2 plaquette and in figure 8 the 3x3 plaquette with and without the matching. 
The values of the bare parameters used are: (3 = 5.662 (pure gauge), (3 = 5.975 - 
K = 0.1676 (n b = 1), and (3 = 6.369 - K = 0.1722 (n b = 2). 

The interesting question is whether such a procedure can be applied success- 
fully to more physical quantities. We have explored the hadron spectrum on a 16 3 • 32 
lattice. In this case the two independent matching quantities have been chosen to 
be the pion and rho masses in order to fix the renormalized quantities at a scale 
appropriate to the problem. The ratio of the two essentially fixes the renormalized 
quark mass and the value of the rho mass, which shows only a mild dependence 
upon the quark mass, the renormalized coupling constant. In reality, they repre- 
sent some complicated definition of the two renormalized parameters. The values 
of the bare couplings needed to obtain a matching of pion and rho masses do not 
differ dramatically from those obtained with the improved perturbative matching 
discussed in the last section: the latter are often the starting point of the matching 
search. We have also found that the correlations of the operators used to extract the 
pion and rho masses at large distances are, at small distances, very useful monitors 
of the matching. These correlations can be obtained in the standard way by invert- 
ing quark propagators and assembling them in the correlation and also extracted 
directly from bermion fields correlations. We have extracted hadron masses using 
pointlike sources with standard and the smeared quark propagators described in 
ref. [7]. 

The results for the hadron spectrum are given in figure 9, 10 and 11 where we plot 
the proton over rho ratio for values of i?2, the square of the ratio of pion over rho 
masses, of ~ 0.68, 0.61 and of 0.50 and of the rho mass in lattice units of ~ 0.705, 
0.70 and of 0.66 respectively. The data points referring to the unquenched case are 
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K = 0.12 


K = 0.16 


quenched 
unquenched 

1 loop 

2 loops 

3 loops 


0.4271(1) 
0.433(1) 
0.4326(2) 
0.4330(4) 


0.4271(1) 
0.4607(5) 
0.4510(1) 
0.4534(2) 
0.4542(7) 



Table 1: Average plaquette for SU(2) with Wilson fermions at (3 = 1.75 and 
K = 0.12 (on a 4 4 lattice), K = 0.16 (on a 8 4 lattice) for the quenched case, 
the unquenched one with rtj = 2 (from HMCA simulations) and the approximated 
results corrected to a given order in the fermion loop expansion. 

taken from ref. [8]. For the lowest value of R2 the unquenched data with approx- 
imately matched values of rho and pion are not available. For the heaviest quark 
mass we show the two unquenched data with the closest matching. In table 2 we 
report the values of the bare parameters used for the matching of theories with dif- 
ferent flavour numbers and those of the corresponding pion, rho and proton masses. 
As one can see the matching is not perfect. In general the value of the proton over 
rho ratio increases when the rho mass or the value of R2 increase: by taking into 
account this effect those points which are not well aligned with the others would do 
it better with a better matching. 

Within our statistical errors and the systematic error due to an imperfect 
matching, there are no sign of physical dynamical fermions effects in the proton 
over rho ratio down to R2 of order 0.5. This conclusion is fully consistent, and 
actually better supported by the smaller statistical errors of our data, with the 
published unquenched data. Recent analysis of accurate quenched data [9] show 
that, by extrapolating to zero lattice spacing and infinite volume, one obtains a 
value for the proton over rho mass ratio close to the experimental value within at 
most ten percent, leaving very little room for dynamical fermion effects. 

In table 3 we compare the values of (3 used for the matching with those that one 
would get for different flavour numbers according to the perturbative eq.(l) where 
the input is the (3 value of the pure gauge case and the output the (3 values of the 
bermion theories. We also show the values of the non integer flavour number which 
would bring the perturbative shift into agreement with the one determined by pion 
and rho matching. In the last two columns of the same table we repeat the same 
exercise by replacing the bare couplings in eq.(l) with those of the Parisi scheme 
given by (3p = f3(P). In this case and for light quark masses the effective flavour 
number agrees with the expected integer value. We notice that the effective flavour 
number corresponding to two bermions {rif = —4) is always roughly twice the one 
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n f 





k 


^conf 


m w 


m p 


M N 


9 / 9 


M N /m p 


-4 


6.4 


0.155 


26 


0.577(4) 


0.705(6) 


1.12(3) 


0.670(10) 


1.59(4) 


-2 


6.1 


0.1557 


20 


0.586(3) 


0.706(7) 


1.120(20) 


0.689(11) 


1.586(25) 





5.767 


0.1582 


20 


0.591(5) 


0.710(7) 


1.125(18) 


0.693(11) 


1.585(24) 


2 


5.5 f 


0.158 




0.568(5) 


0.672(9) 


1.095(20) 


0.714 


1.63(4) 


2 


5.6 


0.156 




0.457(5) 


0.55(1) 


0.87(2) 


0.690 


1.58(7) 


-4 


6.4 


0.158 


36 


0.549(3) 


0.722(7) 


1.144(15) 


0.574(15) 


1.584(22) 


-2 


6.1 


0.158 


40 


0.534(3) 


0.687(8) 


1.079(14) 


0.604(8) 


1.571(20) 





5.7 * 


0.163 


219 


0.562(2) 


0.719(4) 


1.138(10) 


0.611 


1.583(15) 


2 


5.4 f 


0.162 




0.570(10) 


0.717(14) 


1.15(2) 


0.631 


1.60(6) 


-4 


6.463 


0.157 


45 


0.478(4) 


0.663(7) 


1.049(20) 


0.520(10) 


1.582(27) 


-2 


6.1 


0.161 


65 


0.464(4) 


0.658(8) 


1.037(15) 


0.497(7) 


1.576(20) 





5.7 


0.165 


73 


0.460(4) 


0.658(11) 


1.028(18) 


0.489(10) 


1.562(31) 



Table 2: A summary of hadron masses on a 16 3 • 32 lattice used in the extrapolation 
procedure. *Quenched results for (3 = 5.7 and k = 0.163 are obtained by interpolat- 
ing the results for k = 0.160 and k = 0.165 at the same (3 from ref.[8]. ^The results 
for unquenched QCD with 2 dynamical fermions are from ref.[7]. 



n f 




K 


ftpert 


eff 

n f JJ 


r^pert 


eff 
n fPar 





5.7 


.165 


5.7 





5.7 





-2 


6.1 


.161 


6.38 


-1.16 


6.11 


-1.98 


-4 


6.46 


.157 


7.04 


-2.23 


6.48 


-3.95 





5.767 


.1582 


5.767 





5.767 





-2 


6.1 


.1557 


6.46 


-.96 


6.25 


-1.52 


-4 


6.4 


.155 


7.13 


-1.84 


6.74 


-2.9 



Table 3: Comparison between the values of the bare parameters needed for the 
matching and those estimated from eq.(l) 
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of the one bermion case supporting the notion of an "effective flavour number". 

The introduction of negative flavours may have important consequences on 
the U(l) problem in QCD [10]: in particular the contribution proportional to the 
anomaly, splitting the masses of flavour singlet and non singlet pseudoscalars, has 
an opposite sign in the bermion theory with respect to the the fermion theory, 
leading to a splitting in the opposite direction, i.e. lowering the singlet mass. At 
very low quark masses this effect may lead to a massless particle while the pion is 
still massive. We have monitored this phenomenon by calculating flavour singlet 
and non singlet correlations of bermion fields: we have seen a splitting between 
the corresponding masses leading to a lower singlet mass. The simulations with 
almost massless quarks may be problematic due to this effect and may represent 
an ultimate limitation of the method. At the lowest value of the parameter R2 
that we have explored (R2 ~ 0.35) and for n& = 2 where the effect is bigger, we 
have estimated a pion mass of ~ 0.4 and a flavour singlet mass of ~ 0.2 in lattice 
units, away from critical values. All other simulations were performed under safer 
conditions, i.e. at heavier quark masses and/or smaller flavour number. We will 
present elsewhere a detailed analysis of the singlet mass with the bermion method 
as well of other quantities like the light pseudoscalar and vector decay constants 
where a visible contribution of sea quarks is expected to bring the quenched lattice 
data into agreement with experimental values. 

The extrapolations from negative to positive flavour numbers is an economical 
way of estimating dynamical fermion effects. The extrapolation remains smooth 
when the comparison among theories with different flavour content is performed at 
fixed renormalized quantities, i.e. after matching two independent physical quanti- 
ties. 
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Figure 1: The average plaquette for SU(2) at (3 = 1.75 for different values of the 
Wilson parameter K = 0.12,0.14,0.16 as a function of the flavour number. The 
steepest curve refers to the lightest quark mass (K = 0.16). 
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Figure 2: The lxl and 2x2 average plaquettes for SU(S) at (3 = 23.76 and 
K = .1285 corresponding to a negligible (~ 10~ 4 ) quark mass as a function of the 
flavour number. 
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Figure 3: A blow up of the variation of the 2x2 average plaquettes as a function 
of the flavour number, with the same parameters of the previous figure. 
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Figure 4: The variation of the 2x2 average plaquettes as a function of the flavour 
number after the matching procedure described in the text. The values of the lxl 
plaquette and of the perturbative quark mass of ref. [6] are 0.9134 and less than 
10~ 4 respectively. 
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Figure 5: The lxl and 2x2 average plaquettes for SU(3) at j3 = 5.4 and K = 0.162 
as a function of the flavour number. 



15 



0.6 



0.5 



0.4 



0.3 - 



0.2 - 



+ 



+ 



+ 



+ 



0.1 - 



Figure 6: The lxl and 2x2 average plaquettes for SU(S) at approximately constant 
values of the lxl plaquettes and quark mass defined as in ref.[6] of about 0.05 as a 
function of the flavour number, after the matching procedure described in the text. 



16 



0.14 



0.12 



0.1 



0.08 



0.06 - 



0.04 - 



0.02 - 



+■ 



-6 -4 -2 2 4 

U f 

Figure 7: A comparison between the variation with the flavour number of the 2x2 
plaquettes with bare parameters fixed (diamonds) and with the matching (crosses) 
discussed in the text. 
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Figure 8: A comparison between the variation with the flavour number of the 3x3 
plaquettes with bare parameters fixed (diamonds) and with the matching (crosses) 
discussed in the text. 



18 



1.7 



1.65 



1.6 



M^jrrip 1.55 - 



1.5 - 



1.45 - 



1.4 




Figure 9: The ratio nucleon over the p mass as a function of the flavour number at 
approximately constant values of the p and it masses of 0.705 and 0.58 respectively. 
The corresponding numerical values are given in table 1. 
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Figure 10: The ratio nucleon over the p mass as a function of the flavour number at 
approximately constant values of the p and it masses of 0.700 and 0.55 respectively. 
The corresponding numerical values are given in table 1. 
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Figure 11: The ratio nucleon over the p mass as a function of the flavour number at 
approximately constant values of the p and it masses of 0.66 and 0.47 respectively. 
The corresponding numerical values are given in table 1. 
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